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The effect of A site disorder on the magnetic field induced melting of charge order (CO) in half 
doped manganites has been studied using a Monte-Carlo technique. Strong A-site disorder destroys 
CO even without an applied field. At moderate disorder, the zero field CO state survives but has 
several intriguing features in its field response. Our spatially resolved results track the broadening 
of the first order field melting transition due to disorder, explain the unusual dependence of the 
melting scales on bandwidth and disorder, and allow an unified understanding of CO melting across 
all manganites. We also present some results on disorder assisted trapping of met ast able phases in 
the low temperature state as seen in recent experiments. 



I. INTRODUCTION 

The half-doped state in the manganite d"'^ * ^ * ^ * ^ * has been 
thoroughly explored in the recent past. In particular, 
the magnetic field induced melting of the charge order 
(CO) has been probed experimentally to map out the 
hysteretic response, and the bandwidth dependen ce of 
the 'melting' field in systems with lo^^l^, intermediat d^*^ * ^ * 
and strong^ disorder. While these results were primarily 
concerned with the bulk signature of melting, there have 
also been atte mpts to probe the sp atial natur e of field 
melting, direct l}!™™!! and \nd\iect\^^^. 

The main attempts at theoretical understanding have 
been in the 'clean' problem^^^'^^'^ with only an initial 
attemptP^ probing the effect of disorder on field melting. 
The effort to understand the 'clean' problem is under- 
standable because even without disorder there we re many 
issues which have been clarified only recentl} ^^ * ^-*^ * ^ I 
These studies, alongwith an experimental understanding 
of the systematics of disordei^^ set the stage for an 
attempt on the disordered field melting problem. 

To put the experimentally observed disorder effects in 
context, let us recapitulate the physics of the 'clean' sys- 
tem. The bandwidth of a manganite is controlled by the 
mean ionic radius va- At low r^, the half-doped man- 
ganites are insulators with in-plane checkerboard CO, 
dx2_^2 /dy2_^2 orbital order (00), and CE type magnetic 
order. We will simply call this the CE-CO-I phase. The 
large BW materials, with large r^, have a ferromagnetic 
metallic (FM-M) ground state. These two states, and 
an 'A type' antiferromagnetic state, compete with each 
other at X = 0.5. 

For CE-CO-I manganites close to the FM-M bound- 
ary, the application of a modest field converts the ma- 
terial to a FM-M through a first order phase transition 
(FOPT). The critical CO melting field were found to be 
much smaller than the corresponding CO melting tem- 
peratures and is understood to be due to closeness of 
the energies of the competing CE-CO-I and FM-M. The 
melting field increases with decreasing ta as one goes 
further from the FM-M phase. Finally, for very small ta 
systems, the effective single electron bandwidth (BW) is 



so low that the Jahn-Teller (JT) coupling, A say, is able 
to localize electrons by itself and create a CO state with- 
out requiring CE magnetic order. Such charge order is 
naturally impossible to melt with an applied magnetic 
field. The applied field at best would convert the CE- 
CO-I to a FM-CO-I, but not 'melt' the charge order. In 
low disorder systems e.g. Lno.sCao.sMnOa, this broad 
trend of increasing melting fields with decreasing ta has 
been seen^. In our calculation^^, we have shown that 
at intermediate magnetic fields, the ground state ceases 
to be a homogeneous FM-M with density n = 0.5 and 
instead phase separates (PS). The PS is between FM-M 
and AFM-M at larger r^, and between FM-M and FM- 
CO at smaller r^. Further, these coexisting phases are 
all off n = 0.5. Thus, the field melted state is at best a 
percolative metal. 

In real life, while low cation mismatch is possible, as in 
Lno.sCao.sMnOa, eliminating disorder altogether is diffi- 
cult. To this end, there is only one attempt to create 
a disorder free half doped systenP, but it unfortunately 
doesnot probe the nature of the field melted state. In all 
other conventional experiments, the disorder, although 
small is enough to mask the effect above. 

This brings us to a point where we need to carefully 
analyse the effects of weak and strong disorder. First, 
the origin of A-site disorder. In the process of varying 
the r^, one inevitably introduces disorder due to cation 
size mismatch. Disorder makes its appearance in two 
forms, (i) the randomness introduced in the hopping due 
to random deviations of the Mn-O-Mn bond angles and, 
(ii) if the valence state of the dopant is different from 
the parent material, the dopant locations act as scatter- 
ing centers. The extent of structural disorder is quanti- 
fied by the variance, a a, of the distribution of the cation 
radii in Ai-a^A'a^BOsand has been measured for the Ca, 
Sr, and Ba familieJ^^^. From these measurements it 
is known that while in the Ca family (Lno.sCao.sMnOa, 
where Ln=La, Pr, Nd, Sm, Eu, etc.) a a ^ 10~^A^, for 
the Sr family (Lno.sSro.sMnOs) greater size mismatch 
leads to a a ~ 10~^A^. For the Ba family disorder is 
large enough to kill the CE-CO-I state completely even 
without a magnetic fielcP. 
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In our earlier wor M ' ' we classified the 
Lno.sCao.sMnOs family as 'clean' with a disorder 
variance ~ 10~^ A^, while the Lno.sSro.sMnOa 
family was classified as disordered with the disorder 
variance {a a ~ 10~^ A^). We describe our modelling of 
disorder in Section III. In our earlier work we compared 
our results for ideally clean, A = 0, results with the 
Ca family and that of moderate disorder with the Sr 
family. While the A = results do broadly match with 
the trends seen in the Ca family (there is a monotonic 
growth in the critical melting fields with reducing 
va) here we demonstate that A = is not crucial and 
the same trend can be obtained for a window of low A. 
We also track the disorder induced smearing of the first 
order transition. Finally, we show that, starting in from 
A = and going up to strong disorder, the essential 
physics can be rationalized by studying the zero field 
'bandwidth-disorder' phase diagram. 

One is also interested in the spatial nature of the field 
melted state. This line of work has been picked up only 
recently. Experiments on weak and modera tely disor- 
dered system both at half dopin^^^^^^^^^^^^ and away, 
e.g. LPCMCli^J^, show that disorder inhibits the ap- 
proach to equilibrium, leading to non- equilibrium coex- 
istent states with the phase fractions of the coexistent 
states being strongly dependent on the p ath followed to 
bring the system to its final state^S^^fflS]^ present 
work we show results on disorder induced trapping of 
met ast able phase fractions. The detailed nature of the 
nonequilibrium coexistence and its h — T protocol depen- 
dence will be reported elsewhere. 

The paper is organised as follows. In section II we sum- 
marize the key experimental results and in the following 
section III define the model and describe the method for 
solving it. Section IV A. discusses the zero field reference 
state and the effects of disorder on that. Section IV B and 
C discusses the effects of disorder on field melting. and in 
sSection IV D, we present our results on disorder assisted 
trapping of FM-M phase in the low temperature CE-CO- 
I state. Section V concludes the paper. 



II. EXPERIMENTAL RESULTS 

Typically, disorder in manganites ABO3 arise because 
of A or B site substitutions. In this work we consider the 
effects of A-site disordering on the melting of the CO- 
state. We have studied the effects of B-site disorder on 
the CO state at half doping is studied in^^. 

Fig.l.(a)-(b) show the key differences between the 
clean and the disordered results. These plots, recon- 
structed from experimental data^ present the evolution 
of the critical switching fields and CO melting temper- 
atures with decreasing ta (increasing BW) for (a) clean 
(A = 0) and (b) disordered cases. The ± on the critical 
melting fields hco stand for the critical melting field as 
obtained in low temperature magnetic field sweeps in the 
direction of increasing (+) and decreasing (-) field. The 
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FIG. 1: Colour online: The (a) clean ta variation of CO 
melting temperatures (Tco) and magnetic fields (h^o) for 
the Ca family and (b) disordered for the Sr family. Here, the 
± on hco imply forward(-l-) and backward(-) melting fields 
in a low temperature field sweep experiment. The disorder 
variance A in (a) is 10~3 and that is (b) is ~ 1Q~2. These 
plots are reconstructed from experimental datgp^. 



switching fields are defined as the magnetic field value 
where the CO volume fraction switches. This is typi- 
cally concomitant with a sharp change in the resistiv- 
ity (p). The labels on top show the A site dopant and 
(a) is for the Lno.sCao.sMnOs family and (b) is for the 
Lno.sSro.sMnOa family. Increase in A (equivalently de- 
crease in va or reduction in BW) for the 'clean' system 
Figl (a) increases the CO melting temperatures and the 
switching fields h^Q. 

For a range of va 1.36-1.32, while there is a increase 
in Tco and h^Q , the corresponding quantities in the dis- 
ordered case initially rise a bit and then is very strongly 
suppressed and drops to zero at va ^ 1.32, which cor- 
responds to SmcsSrcsMnOgP. For materials with larger 
disorder such as the Ba family, the CO state does not 
exist even at /z = cP. 

While these measurements probe the bu lk effects of 
disorder, recent experiments have also probecP^l^m^n^^ 
the effects of disorder on the spatial nature of the melt- 
ing phenomenon. In particular (Chaddah et.al. in their 
work^^). Fig. 2 (a) show that even in the weak disorder 
samples there is marked non-equilibrium phase coexis- 
tence. Here samples of Lao.sCao.sMnOa are field cooled 
from 320K to 5K at different fields. The sample cooled 
in higher field has greater fraction of FM-M regions co- 
existing with AF-I regions. However, as shown, for the 
6T field cooled sample, a small increase in temperature 
to about 80K, after isothermally reducing the applied 
field to IT (at 5K), makes the system lose this FM vol- 
ume fraction. Infact, the magnetization drops to that 
of the 1 Tesla field cooled sample at the same tempera- 
ture. While magnetization by itself is not a good probe 
of amount of FM regions (which domains can be mis- 
aligned), the corresponding resistivity shown in Fig.2.(b), 
too show an increase in the resistivity for the (6 Tesla) 



3 




100 200 300 100 200 300 
T(K) T(K) 



FIG. 2: Colour online: (a) The variation of magnetization^^ 
with temperature for samples cooled in IT and for sample 
cooled in 6T field, (b) The corresponding resistivities. The 
red line corresponds to the warming path in IT field. The 
multivalued nature of the magnetization and the resistivity is 
seen in the varied values of these quantities when the systems 
are cooled in different fields and then isothermally brought to 
the same (/i, T) point. The non-equilibrium nature of the (6T) 
field cooled system is seen in the fact that the magnetization 
and the resistivity both approach the IT field cooled results on 
increasing temperature to 80K, after decreasing the field (on 
the 6T sample) to IT isothermally at low temperature(5K). 
The thermal hysteresis in the IT case (in b) is due to the 
accompanying metallic to insulator first order transition with 
decreasing temperature. 



Here, c and are annihilation and creation operators for 
eg electrons and a, jS are the two Mn-e^ orbitals dx2_y2 
and dsz^-r"^^ labelled (a) and (b) in what follows, are 
hopping amplitudes between nearest-neighbor sites with 
the symmetry dictated form: t^^ = f^^ = t, t^^ = t^^ = 
^3, = t^^ = -t/V3, tl^ = tl^ = t/V3, where x and y 
are spatial directions We consider effectively a lattice of 
Mn ions and treat the alloy disorder due to cationic sub- 
stitution as a random potential at the Mn site picked 
from the distribution PAi^i) = ^(^(e^ - A) + S{ei + A)). 

The eg electron spin is af = 4aa^aa'^ic,^' ^ where 

the F's are Pauli matrices. It is coupled to the t2g spin 
via the Hund's coupling J//, and we assume Jn/t':^ I. A 
is the coupling between the JT distortion = {Qix^ Qiz) 

and the orbital pseudospin r^^ = Y^^^ AacT^^(3^il3cr^ 
K is the lattice stiffness. We set t = 1, = 1, and treat 
the Qi and as classical variables. The chemical po- 
tential ji is adjusted so that the electron density remains 
n = 1/2 which is also x = 1 — n = 1/2. Note that this 
modelling of disorder as random fluctuating potentials is 
based on the fact the while hopping disorder is crucial to 
the formation of spin-glass phases, it is onsite disorder 
rather than hopping disorder that plays the dominant 
role in the kind of physics we are interested irPl. 

B. Method 



field cooled sample at about 80K and essentially coincid- 
ing with that of the 1 Tesla field cooled sample at this 
temperature. These results imply the non-equilibrium 
nature of the low T state and the tunability of phase 
volumes depending on h — T protocol. 

We will discuss our results in context of the trapping 
of the FM-M phase in the (CE-CO-I) ground state and 
its systematics with changing disorder, while we leave 
the issues of the nonequilibrium nature and the protocol 
dependence for a later work. 

III. MODEL AND METHOD 
A. Model 

We consider a two band model for eg electrons, Hunds 
coupled to t2g derived core spins, in a two dimensional 
square lattice. The electrons are also coupled to Jahn- 
Teller phonons, while the core spins have an AF superex- 
change coupling between them. These ingredients are all 
necessary to obtain a CE-CO-00 phase. We include the 
effect of disorder through an on site potential. 

{ij)a i i 

^Jaf ^^'^j -^Y. ^^ ^^ + y - ^ 5Z ^ 

{ij) i i i 



We employ a real space exact diagonalization (ED) 
based Monte Carlo (MC) techniqu^^ that allows us to 
perfrom calculations on system sizes of up to 40^. Details 
are given Typically results are averaged over 20-30 
disorder realizations. 



C. Physical quantities 

In order to study the evolution of the system with 
applied magnetic field we track various physical quan- 
tities in real space and momentum space. We com- 
pute the 'one point' distribution of lattice distortions, 
P{Q) = Hi^{Q - Qi)^ where Qi = |Qi|, spatial Qi.Qj 
correlations, DQ{q) = 'Yij{QiQj)e'^^'^^^~^^\ and spin- 
spin correlations, S{q) = - Angular 
brackets represent a thermal average. We also compute 
the volume fraction of the charge ordered region in the 
lattice from direct spatial snapshots of the charge distri- 
bution. To measure the volume fraction, we tag a site 
with a particular color if the site has n > 0.5 and is sur- 
rounded by the four nearest neighbor sites with n < 0.5 
and vice verse (i.e. a site with local anti-ferro-charge 
correlation is marked with a particular color). Simi- 
larly, if the difference between the charge density at a 
site with its nearest neighbours is less than a threshold, 
that site is tagged by a different color, ie., the charge 
^uniform regions are marked by this color. For inter- 
mediate cases, we use an interpolative colour scheme. 
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The volume fraction is necessary for studying inhomo- 
geneous melting where the momentum space structure 
factors are not a good measure of the amount of CO 
in the system. Further, the spatial snapshots of the real 
space charge density also directly provide visual informa- 
tion on the melting process. While the indicators above 
measure the spatial correlations and spatial evolution, 
the metallic or insulating character is tracked via (low 
frequency) conductivity, (Jdci and the density of states 
(DOS), N{uo) = (-^ 5{ui) — en)), where are the elec- 
tronic eigenvalues in some MC background and the an- 
gular brackets indicate thermal average. We track all 
the above quantities as a function of temperature and 
applied magnetic fields for studying the CO melting phe- 
nomenon. Results are disorder averaged over 15-20 dis- 
tinct disorder realizations. 



IV. RESULTS 

We begin by looking at the effects of disorder without 
any applied field. We will then look at the combined 
effects of both disorder and magnetic fields and discuss 
the fate of the charge-orbital-spin ordered state. 



A. Zero field 

We first look at the zero field disorder effects on the 
charge ordered state by cooling the system in zero field at 
different disorder strengths for various values of electron- 
phonon coupling(A). Our goal here is to probe the effects 
of disorder on the charge ordered state using bulk indica- 
tors like the volume of the CO state (Vco) indicating the 
extent of the CO in the ground state, the charge order 
structure factor, giving the coherence of the CO regions 
and analysis of spatial snapshots of the charge density 
fields {ui). 

Figure 3 (a) and (b) show the disorder averaged volume 
fraction of the charge ordered regions in the ground state 
as a function of disorder strengths for various A values 
indicated in (a) and the corresponding tt, tt charge order 
structure factors in (b). Both for small A ~ 1.45 and for 
large A ~ 2.0, the volume fractions of the charge ordered 
regions decrease rapidly with A, the disorder strength, 
while for intermediate A ~ 1.6 — 1.7, the volume fraction 
of the CO regions remains robust till A ~ 0.12, with the 
A = 1.7 case having its CO volume fraction suppressed 
slightly more rapidly than that for A ~ 1.6 case. The 
CO structure factor provides a clearer picture of the way 
disorder destabilizes the CO state. For weak and strong 
A cases till about A = 0.1, while the volume fraction of 
the CO regions is about 80%, the corresponding structure 
factors diminish rapidly, indicating formation of domain 
walls between different CO regions. For the intermediate 
A case, the system presents a single phase CO state till 
A ~ 0.12, beyond which the A = 1.7 case is more prone 
to losing the CO state than the A = 1.6 case. This fact is 




FIG. 3: Colour online: The variation of (a) volume fraction 
of charge order(CO) and (b) the long range nature of the 
CO state, i.e., the CO structure factor Z)g(7r,7r), with disorder 
strength at various electron phonon coupling (A). Note the 
fall in the Vco and Dq{7v,7v) with disorder is most rapid for 
the low and large A cases. Also note that long rangeness 
of the CO regions is lost much more rapidly with increase 
in disorder than the volume fraction of the CO itself. All 
results are obtained for 16^ systems and are averaged over 
many disorder realisations. 



seen exemplified in the spatial snapshot analysis to which 
we turn next. Figure 4 shows the snapshots of the CO 
state {rii) at three values of disorder, from weak, through 
intermediate, to strong from top to bottom, for the four A 
values increasing from left to right (see the corresponding 
caption for the values). Figure 5 shows the snapshots of 
the corresponding magnetic states. 

(i) Small and large A regimes: At small disorder, the 
low A case (the extreme left column) begins developing 
metallic regions, which are concomitant with local dis- 
ruption of the CE chains. With increase in disorder, these 
metallic regions grow along with appearance of line-like 
AF states. Further increase of disorder to ~ 0.15 leads 
to formation of distinct FM-M regions. The strong cou- 
pling case, which essentially consists of almost site local- 
ized electrons with weak overlap with neighbouring sites, 
starts forming clear domain walls at low and intermedi- 
ate disorder and finally at large disorder, forms a badly 
disordered polaronic state. The corresponding magnetic 
state is intially disrupted by formation of local CE regions 
of opposite handedness, which at large disorder devel- 
ops patches with local (tt, tt) magnetic correlations. The 
small and larger A states under the effect of disorder can 
be understood in terms of disorder hindering the growth 
of the ideal CE-CO-I state from the intermediate temper- 
ature states. For low A, at intermediate temperature and 
A = 0, the system is in a FM-M phase which in absence 
of disorder gives way to a CE-CO-I state at low temper- 
ature. However in the presence of disorder, a fraction of 
the FM-M phase remains trapped at low temperatures. 
Similarly, for large A case, since the effective charge or- 
der stiffness is small at large A values^^, the intermediate 
state is a disordered polaronic phase and it is this state 
which remains partially trapped due to disorder when the 
system is cooled to low temperatures. 





FIG. 4: Colour online: The spatial variation of the electronic 
charge densities for four A values (1.45,1.6,1.7,2.0), from left 
to right, at weak (0.08), intermediate (0.10) and strong(0.15) 
disorder strengths, from top to bottom. 



Since the metallic regions at low T can gain kinetic en- 
ergy by delocalizing electrons within the regions, double 
exchange promotes FM order and the system develops 
FM-M patches, for the low A case. For large A case how- 
ever, the since the on site localization tendency (which 
grows as A^) is larger than the gain from kinetic energy, 
the double exchange is suppressed and antiferromagnetic- 
superexchange dictates a local (7r,7r) antiferromagnetic 
regions existing in a matrix of small CE regions (which 
exists in regions with local CO). 

(ii) Intermediate A regime: Increasing A from weak to 
intermediate values lead to a deepening of the CE-00- 
CO minima in the energy landscape small, where, disor- 
der is unable to trapthe system partly in the FM-M phase 
(which in a Ginzburg-Landau scenario would no longer 
remain metastable at low temperatures). Similarly de- 
creasing A from strong couplings increases the effective 
CO stiffness, thus pushing up, in energy, the disordered 
polaronic phase, which again for weak and intermediate 
disorder would be unstable at low temperature. 

The trend in stability of the CO state with A suggests 
a possible intermediate A where the CE-OO-CO state is 
most robust to disorder. This corresponds to A ~ 1.6, for 
the parameters used in the present work. Results at A = 
1.7 are specifically shown to demonstrate the weakening 
of the CO state compared to A = 1.6. 

(iii) A — A phase diagram: Fig. 6. (a) shows the variation 
of Too with A at various A. 6.(b) shows the A — A phase 
diagram. In Fig. 6. (a), the low disorder Too increases 
with increasing A, but for A ~ 0.1, the Tco at larger 
A starts falling and by A ~ 0.12, goes below those for 
weaker A values. 

To understand this behaviour, we mapped out the the 
A — A phase diagram shown in Fig.6.(b). Since at very 
low A and A = 0, the system is in a FM-M state, at small 
disorder, this state will be a disordered FM-M. Thus, the 



FIG. 5: Colour online: The real space snapshots of the mag- 
netic state corresponding to the charge density profile shown 
in Fig. 4. The parameter points are the same as in Fig. 
4. (See the caption of Fig. 4). Minus one(green) indicates 
perfect antiferromagnetic bonds and plus one (red) indicates 
perfect ferromagnetic bonds. 



CO region in the A — A parameter space, on the lower 
A side is bounded by a FM-M at zero disorder and by 
a disordered FM-M at intermediate disorder. On the 
other hand, at large A and zero disorder, the system, at 
T=0 will remain in the CO state at all finite A values 
(ignoring quantum fluctuations). However, at any finite 
T, small thermal fluctuations will destroy the long range 
nature of the charge order, although the electrons will 
be almost site localized. This will happen due to the 
essential on site localization of the electrons which leads 
to suppression of the effective charge order stiffness. By 
the same token, at any finite A, there is always a critical 
A, beyond which the system, even at T=0, turns into a 
disordered polaronic state. 

The asymptotic limit where A is large enough to allow 
perturbation in kinetic energy can be analytically han- 
dled. Infact this very general disorder induced suppres- 
sion of long range order is well understood in classical 
models. For the random field Ising model (RFIM), this 
has been studied in detail^^. The real regime where the 
materials exist, A ~ 1.5 — 1.7, is non-perturbative and has 
to be understood as the regime interpolating between the 
large A asymptote and the small A < 1.4 FM-M. This nat- 
urally causes, as seen in Fig. 6(b), the region for stable 
CO to be finite at any finite value of A and explains why 
with increasing A, the Tco falls more rapidly at finite 
disorder. Further, the two kinds of CO regions, CD-I 
and CD-II, in the Fig. 6(b), refer to disordered FM-M 
and disordered polaronic regime with AF magnetic back- 
ground. The dashed lines on Fig. 6 (b) represents param- 
eter points along which calculations were done. We shall 
later use the A — A phase diagram to explain experimen- 
tal observations on the material systematics of the critical 
melting fields. 
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FIG. 6: Colour online: (a) The variation of Tco with disorder 
strength A, at various electron-phonon couplings, (b) The 
A — A phase diagram. Note, the charge disordered regions 
at weak A < 1.45 is very different from that which occurs 
at larger A > 2. . The weaker A charge disordered phase 
has uniform charge density and the system is a ferromagnetic 
metal. At large A, the system is essentially site localized and 
the nearest neighbour overlap is negligible and the system can 
lose the charge order under small perturbation and turn into 
a disordered polaronic state. 



B. Disorder induced rounding of first order 
transition 
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FIG. 7: Colour online: The comparison of various physical 
quantities Vco, Sfm and p for the clean and the disordered 
cases. The sharp first order like transitions, in (a), in the 
clean case are replaced by a smoother and 'rounded' transition 
in Vco (blue curves) in (b). The accompanying Sfm ferro 
magnetic order parameter (red curves) reaches its maximum 
(unity) through a series of small steps in (b) as opposed to a 
sharp jump in (a). The resistivity (black curves) in (d) also 
similar behaviour of changing in small steps in contrast to the 
sharp change in (c). 



When a first order phase transition (FOPT) occurs in 
a disordered background, the transition ceases to be dis- 
continuous and acquires a smooth character. Here we 
track the magnetic field induced FOPT between the CE- 
CO-I and the FM-M as it takes place in the disordered 
background and compare to what we found in the clean 
case^^. Fig. 7 and 8 show this evolution of the CO state 
as we sweep the magnetic field after cooling the system 
to low temperature. While Fig. 7 tracks measures of 
bulk properties Vco^ Sfm and p, as a function of applied 
field in the clean (a), (c) and disordered (b) ,(d) cases, 
Fig. 8 looks at the actual spatial profile of the CO melt- 
ing. The spatial profile in Fig. 8, for the clean system 
(top panel) is seen to resist the field up to h/t = 0.06, 
beyond which it abruptly goes to a percolative FM-M 
state. The disordered case (bottom panel) shows a more 
gradual trend in the melting, it starts by creating small 
metallic regions which grow over a window of field values 
to reach the final percolative metallic state. This gradual 
loss of CO volume fraction is corroborated by Fig. 7(b) 
and the abruptness of the transition on the clean case is 
correlated with Fig. 7(a). The fact the the systems are 
metallic is seen from the resistivities shown in Fig. 7(c) 
and (d) which follow the abruptness of the change in Vco 
in the clean case and smoothened change in disordered 
case. 



C. Field melting & the relation to the A — A phase 
diagram 

Here we look at the detailed behaviour of the ther- 
mal and magnetic melting scales for the CO state with A 
and disorder. Since it would also help making compari- 
son with experiments, we restate the broad experimental 
results. While for the systems with cta ^ 10~^A^ (e.g. 
the Ca family), the magnetic melting scales increase with 
decreasing r^, for systems with a a ~ 10~^A^ (e.g. the 
Sr family), the melting scales initially increases with de- 
creasing ta and then are strongly suppressed and even- 
tually go to zero with further decrease in r^. These are 
shown in Fig. 1(a) and l.(b). 

We follow the same h^T protocols as in the experi- 
ments, ie, we sweep up and down in h at low T after 
cooling at = 0. In terms of the numerics we start from 
zero field and then first increase the magnetic field up 
to h/t ~ 0.2 and then reduce to zero in steps of 0.01. 
As discussed earliei^^, the forward melting field is de- 
fined as h^Q and that in the downward sweep is h'^Q. 
Fig.9(a)-(d) show the results of our calculation starting 
with A/t = in (a) to strong disorder A ~ 0.12 in (d). 
We clearly see that apart from the expected overall sup- 
pression in melting scales with disorder, there is a gradual 



FIG. 8: Colour online: The spatial snaphots of charge densities for zero field cooled samples, which are then subjected to 
a sweep in the magnetic field at low temperature for clean (Top panel) and disordered ( A /t = 0.12) (Bottom panel). The 
snapshots are shown for field values (h/t=0. 04, 0.05, 0.06, 0.07, 0.08, 0.09) from left to right, for the increasing part of the field 
cycle. The corresponding CO volume fractions are shown in Fig. 7(a) and (b). Clearly, in the clean case, the CO remains stable 
up to h/t ^ 0.06 and then collapses within a window of A{h/t) = 0.01 to a state with about 30% CO regions and 70% FM-M. 
The disordered case however, starts losing CO dX h/t — 0.05 and collapses, within a window of A(h/t) = 0.04 into a percolative 
metallic state with roughly similar composition as in the clean case. 



downturn with increasing A (or decreasing va)- In (a) the 
melting field diverges, i.e., the CO state is stable without 
CE order. In (b) the melting scales become finite in the 
same A regime and begin to get strongly suppressed in 
(c) and (d). 



To explain the results we refer back to Fig. 6(b), the 
A — A phase diagram. Here we have marked by the red 
dashed lines the (A, A) combinations for which we have 
shown results in Fig. 9(a)-(d) (with progressive increase in 
the values of disorder from a to d). From Fig 6(b) it is 
clear that while at A = 0, the system has only one bound- 
ary with the CD-I region, so the melting scales are the 
smallest at A ~ 1.45 and increase with increasing A. For 
all other dashed lines the systems encounter boundaries 
with CD-I and CD-2, so there is a generic suppression 
in the magnetic and thermal melting scales. Thus the 
presence of disorder causes a competition between the 
long range order of the CO state and random pinning ef- 
fects of disorder which weakens the effecting CO stiffness 
and hence the melting temperatures and melting fields^ ^. 
This at large A and presence of the FM-M at low A places 
the CE-CO-I in a small Awindow. While the Ca family 
wold fall in the low A category, the Sr family would, 
in the same A regime, be placed at A ~ 0.12. The Ba 
family would be beyond A ~ 0.15. The A — A phase 
diagram gives a natural framework to rationalize both 
the clean and disordered field melting systematics in an 
unified manner. 



D. Disorder induced trapping of metastable states 

We now briefly turn to certain aspects of disorder in- 
duced coexistence in the half doped manganites The ex- 
perimental results as discussed in Section II have three 
major aspects, one the disorder induced trapping of 
metastable phase in the low temperature state of the sys- 
tem, two, the non-equilibrium nature of such coexistent 
states and three, the strong protocol dependence of the 
phase fractions. Here we present our results only on the 
first aspect, while we will report the remaining two in a 
later work. 

The results shown in Fig. 10 are obtained by cooling 
the systems in presence of various fixed magnetic fields 
at different disorder strengths. All results are on 16^ 
systems and the results are averaged over many disorder 
realizations. We use two indicators for the study, one, is 
the volume faction of the CO regions (Vco) and the other 
is the measure of excess number of FM bonds in the sys- 
tem over that which exists in the CE phase (AVfm)- In 
the CE phase, each site in the plane has four bonds two 
of which are ferromagnetic and two antiferromagnetic. 

Thus if we define AVfm = ^""Ifm^ , then in the CE 

CE 

phase AVfm = 0, in the FM phase it is 1 and is negative 
if the number of bonds antiferromagnetic bonds are more 
than the number of ferromagnetic bonds. In the case of 
G-type AF phase it is -1. Although AVfm cannot distin- 
guish between two spatial patterns with the same number 
of AF and FM bonds, it is a useful measure of net FM 
volume on the system and in the present definition gives 
the net increase or decrease in the FM volume fraction 
from that which exists in the CE phase. 
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FIG. 9: Colour online: (a) to (d) show the variation of the disorder averaged Tco, and h'^ for A = 0,0.06,0.1,0.12 

respectively. While for A = 0, the melting fields diverge beyond a critical A and the Tco grows, in the shown A window, 
increasing the strength of disorder A not only suppresses the overall melting scales, but also changes the A (or ta) dependence 
qualitatively. Within the same A window, increasing disorder causes a gradual downturn in all the melting scales at weak 
disorder and finally strongly suppresses these scales at strong disorder A ^ 0.12. Error bars for the strong disorder case (d) 
are the largest and are given as an estimate of the maximum error in the numerics. 
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FIG. 10: Colour online: (a) The evolution of CO volume 
fraction with h for varying A. The finite Vco at /i ~ 0.05 
and beyond is due to PS at intermediate magnetic fields. For 
lower A the system loses most of Vco at low fields before 
PS into FM-M, FM-CO and AF-M. For A - 1.6, the system 
directly PS between AF-M and FM-CO at /i - 0.05, without 
the occurrence of the low field FM-M state, (b) The effect of 
disorder on Vco- Results are for A = 1.5 and A = 0.05, 0.16. 
When compared to A = (dashed line), increasing disorder, 
at low /i, leads to larger Vco- At larger /i, however, this 
effect is masked by the field induced PS. (c) The weak field 
response of the the magnetic ground state at A = 0. For 
h < 0.012, the system remains in the CE state. Beyond this 
within a Ah ^ 0.002, the system loses about 10% of the 
antiferromagnetic bonds to ferromagnetic ones, (see text for 
a definition of AVfm)- (d) Remaining at h < 0.012, where 
the clean ground state is CE-CO-I, this shows the trapping of 
the higher temperature FM with increasing disorder. Results 
are shown for three cooling fields. 



(i) Effects of magnetic field in a clean systems: We 
begin by looking at the effect of magnetic field on Vco 
when there is no disorder in the system. Fig. 10(a) 
shows the variation of Vco with applied magnetic fields, 
at three A values. It was shown in^^, that the system 
at half doping phase separates, at finite magnetic fields, 
into 'FM-M AF-M' at low A and between 'FM-CO AF- 
M' at A > 1.6. Here we see that for A < 1.6 the system 
at low fields {h ~ 0.02 — 0.04), actually lose almost all 
the CO regions to a FM-M regions and then regain a 
certain volume fraction of the CO state at intermediate 
field values. For A > 1.6, the system directly goes into a 
phase separated state, without the low field FM-M phase. 
Fig. 10(c) shows AVfm^ for low fields, for A = 1.5. For 
h ~ — 0.01, the system resists any significant change 
in the CE- phase. At h = 0.012, it converts about 10% 
of the AF bonds to FM bonds. Further increase of field 
makes the system lose more AF bonds. However, we will 
restrict our attention to h < 0.01 and look at the effects 
of disorder on the CE phase. 

(ii) Interplay of magnetic field and disorder: Fig. 10(b) 
shows Vco as a function of magnetic field at two disor- 
der values for a typical low A(= 1.5) case. While beyond 
h = 0.05, the phase separation tendency masks the dis- 
order effects, we focus our attention to low field. For 
h < 0.05, we see that while for those field values where 
the clean system is in the CO state, increasing A reduces 
the volume fraction of the CO regions and for those values 
of h where the clean system would be a FM-M, disorder 
traps some fraction of the system in the CO state. This 
clearly shows the generic effect of disorder assisted trap- 
ping of the nearby (in energy) metastable phases in the 
low temperature state of the system. Fig. 10(d) shows 
the variation of AVfm with disorder at various magnetic 
field values. At very low fields h = 0.002, the system 
presents no change to the CE phase when there is no dis- 
order, however, in the presence of disorder, the ground 
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state has a certain volume fraction of FM region and this 
fraction grow to about 5% as the disorder strength A 
grows from to 0.2. At larger A ~ 0.3 the system does- 
not have any clear order due to strong random pinning 
effects of disorder. Further, increasing magnetic field and 
disorder allows one to achieve a continuum of values for 
AVfm^ which together with the results in Fig. 10(b), 
shows the tunability of the CE-CO and FM-M volume 
fractions by cooling in different applied field^. 

V. CONCLUSIONS 

Using a real space Monte-Carlo scheme we have shown 
how the charge ordered state in the half-doped mangan- 



ites is affected by the presence of disorder and an ap- 
plied magnetic field. We illustrated the disorder induced 
broadening of the melting transition and mapped out 
the full h — T phase diagram for various combinations 
of bandwidth and disorder. We could explain the coun- 
terintuitive bandwidth dependence seen in the Sr based 
manganites. We motivated the trends seen in the detailed 
numerical results by analysing the stability of the h = 
charge ordered state in response to disorder. Finally, 
we observe that disorder leads to trapping of metastable 
phases in the low temperature state. Different field cool- 
ing pathways will trap different amounts of FM-M in the 
nominal 'CE-CO-I' ground state. 

We acknowledge use of the Beowulf cluster at HRI. 
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